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Nitroglycerine  has  been  widely  used  as  an  ingredient  of  explosives  and  propellants  for  a  long 
time.  It  is  produced  by  the  glycerine  nitration  reaction.  Many  fire  and  explosion  incidents 
have  occurred  during  its  manufacturing  process  owing  to  its  severe  thermal  instability  in 
both  the  nitration  reaction  and  the  purification  processes.  The  kinetic  parameters  of  this 
reaction  and  the  criteria  for  stable  production  have  not  so  far  been  clearly  identified.  In  this 
investigation,  we  estimate  its  reaction  kinetic  parameters  from  the  experimental  results 
obtained  using  a  tubular  flow  reactor  by  a  numerical  analysis  method.  These  reaction  kinetic 
parameters  are  then  used  to  evaluate  the  safe  operating  conditions  in  a  CSTR  for  the  Biazzi 
continuous  nitroglycerine  production.  The  results  are  important  and  useful  in  the  produc¬ 
tion  process  of  nitroglycerine. 

©  2007  The  Institution  of  Chemical  Engineers.  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nitroglycerine  (NG),  also  known  as  glycerol  trinitrate,  is 
usually  produced  from  the  glycerine  nitration.  Its  nomencla¬ 
ture  is  1,2,3-propanetriol  trinitrate,  chemical  formula  is 
C3H5(0N02)3,  and  molecular  weight  is  227  g/mol.  The  indus¬ 
trial  product  is  a  yellowish  oily  liquid  having  a  density  of 
1.591  g/cm3  at  25  °C.  Nitroglycerine  has  explosive  properties 
and  has  been  widely  used  as  ingredient  of  explosives  and 
propellant  for  more  than  a  hundred  years.  It  is  the  main 
component  of  high  explosives  such  as  dynamites,  as  well  as  an 
ingredient  in  most  mining  explosives.  It  was  used  extensively 
in  smokeless  powder  and  together  with  nitrocellulose  as  a  so- 
called  double  base  propellant,  or  with  nitrocellulose  and 
nitroguanidine  as  a  triple  base  propellant,  which  are  used  in 
larger  calibre  projectiles.  Furthermore,  it  has  been  used  as 
ingredient  of  the  solid  propellant  in  rockets. 

Nitroglycerine  has  very  high  impact  sensitivity  and  thermal 
instability.  Thus,  many  accidents  have  occurred  during  its 
manufacture  and  transport.  According  to  Biasutti  (1985),  500 
major  accidents  occurred  in  the  explosive  industry  worldwide 
during  the  period  from  1769  to  1980.  Among  these  accidents,  121 
involved  nitroglycerine  which  is  approximately  24.2%  of  the 
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total.  Han  et  al.  (1988a)  show  that  there  were  43  industrial 
accidents  involving  nitroglycerine  in  the  period  from  1954  to 
1982.  Obviously,  the  manufacturing  process  of  nitroglycerine  is 
more  dangerous  than  other  chemical  processes. 

The  manufacturing  technique  for  nitroglycerine  includes 
both  batch  and  continuous  processes.  Typical  examples  of  the 
batch  method  are  the  Nobel  process  from  1862  and  the  Nathan, 
Thomson  and  Rintoul  process  from  1908.  A  continuous  method 
was  first  developed  by  the  Schmid  process  in  1927  followed  by 
the  Biazzi  process  in  1935.  Subsequently  Nilssen  and  Brunnberg 
developed  a  new  injection  nitration  process  in  1950.  In  1978  the 
Hercules  Company  in  U.S.A.  designed  another  new  tubular 
process  (Urbanski,  1965;  Kaye,  1978a). 

Although  the  manufacture  of  nitroglycerine  has  been 
developed  for  over  140  years,  the  dangerous  industrial 
procedures  involved  have  almost  remained  the  same.  These 
procedures  include  preparing  of  mixed  acid,  nitration  of 
glycerine,  separation  of  product  and  wastes,  washing  of 
product,  filtration  of  waste,  recycling  of  nitroglycerine  from 
waste  acid  and  treating  of  waste  acid.  Most  of  industrial 
accidents  have  occurred  during  the  glycerine  nitration  stage  or 
the  separation  stage  (approximately  50.4%,  Biasutti,  1985). 
Therefore,  improvements  in  manufacturing  procedures  and 
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Nomenclature 

A 

pre-exponential  factor  of  the  Arrhenius  equa¬ 
tion 

C 

concentration 

CP 

molar  heat  capacity  of  individual  component 
(kj  mor1  FT1) 

Cv 

heat  capacity  of  reactant  (kj  l-1  K-1) 

ACpd 

specific  heat  of  mixed  acid  at  standard  state 
(kj  mor1  FT1) 

E 

activation  energy  (kj  mol-1) 

h 

overall  heat  transfer  coefficient  of  ambient 
cooling  medium  (kj  m-2  K-1  min-1) 

H 

enthalpy  (kj  mol-1) 

AHd 

heat  of  dilution,  kj  mol-1 

A  Hr 

heat  of  reaction  (kj  mol-1) 

k 

reaction  rate  constant;  dimension  depends  on 
kinetics 

m 

reaction  order  of  nitric  acid 

M 

initial  concentration  ratio  of  nitric  acid  to  gly¬ 
cerine 

n 

reaction  order  of  glycerine 

Qg 

dimensionless  heat  generation  rate;  definition 
in  Eq.  (61) 

Or 

dimensionless  heat  removed  rate;  definition  in 
Eq.  (62) 

r 

reaction  rate;  dimension  depends  on  kinetics 

R 

universal  constant  (kj  kg-1  mol-1  K-1) 

RF 

weight  ratio  of  mixed  acid  to  glycerine 

Rn 

percentage  of  nitric  acid  in  anhydrous  mixed 
acid 

Rw 

percentage  of  water  in  mixed  acid 

s 

the  external  surface  area  of  cooling  coils  (m2) 

T 

temperature  (K) 

Vo 

flow  rate  of  reacted  solution 

V 

volume  of  reactor  (1) 

wm0 

inflow  weight  of  mixed  acid  in  reactor  (g) 

Wmo 

inflow  weight  rate  of  mixed  acid  in  reactor 
(kg  mirr1) 

Wmf 

outflow  weight  of  mixed  acid  in  reactor  (g) 

Wmf 

outflow  weight  rate  of  mixed  acid  in  reactor 
(kg  min-1) 

X 

fractional  conversion  of  component 

Greek  symbols 

P 

dimensionless  ambient  temperature  in  cooling 
system;  definition  in  Eq.  (50) 

8 

definition  in  Eq.  (57) 

00 

definition  in  Eq.  (55) 

01 

definition  in  Eq.  (56) 

<P 

definition  in  Eq.  (52) 

0 

definition  in  Eq.  (51) 

definition  in  Eq.  (53) 

#ref 

definition  in  Eq.  (54) 

r 

resident  time 

Subscripts 

a 

in  coolant  system 

f 

final  condition 

G 

C3H5(OH)3 

j 

index  of  components 

N  HN03 

NG  C3H5(0N02)3 

o  initial  condition 

ref  standard  state 

S  H2S04 

W  H20 


the  development  of  safety  control  techniques  are  important 
requirements  for  this  reaction. 

The  kinetic  parameters  of  the  nitroglycerine  reaction  and 
the  operating  conditions  for  its  industrial  manufacture  are  not 
found  in  the  open  literature,  presumably  because  of  its 
military  importance.  In  this  investigation,  we  have  estimated 
the  kinetic  parameters  of  the  nitroglycerine  reaction  from  the 
experimental  results  obtained  in  a  tubular  flow  reactor  by  a 
numerical  analysis  method.  Furthermore,  we  evaluate  the 
safe  operating  conditions  in  a  CSTR  of  the  Biazzi  continuous 
process  using  these  kinetic  parameters. 


2.  Evaluation  of  kinetic  parameters  of 
glycerine-nitric  acid  reaction  system  from  the 
experimental  results  in  a  tubular  flow  reactor 

The  elemental  reaction  of  producing  nitroglycerine  is  ester¬ 
ification  reaction  of  glycerine  and  nitric  acid.  This  chemical 
reaction  can  be  expressed  as: 

C3H5(OH)3  +  3HN03H^S4C3H5(0N02)3  +  3H20  (1) 

In  this  chemical  reaction,  the  mixed  feed  acid  contains  both 
HN03  and  H2S04.  The  reagent  H2S04  acts  as  a  dehydrating 
agent  (Ren,  1994).  Therefore,  only  the  HN03  participates  in  this 
reaction.  The  chemical  consumption  can  be  expressed  by 

CG  =  CG,o(l-xG)  (2) 

and 

Cn  =  Cn,o  -  3CG;oXG  =  CG;o(M  -  3xg)  (3) 

where  CG  and  CN  denote  the  concentration  of  glycerine  and 
nitric  acid,  respectively.  xG  represents  the  fractional  conver¬ 
sion  of  glycerine  and  M  is  the  initial  concentration  ratio  of 
CNj0  to  CGf0.  Assuming  the  reaction  is  n  order  in  the  glycerine 
and  m  order  in  the  nitric  acid  and  that  the  reaction  rate 
constant  can  be  expressed  as  Arrhenius  form,  where  A  and  E 
denote  the  frequency  factor  and  activation  energy,  respec¬ 
tively.  Then,  the  reaction  rate  of  glycerine  and  nitric  acid  can 
be  expressed  as 

-rG  =  kCSCS=Ae-E/RTQCS 

=  A  e-E/RTQ+m(l  -  xG)n(M  -  3xG)m  (4) 

2.1.  Mass  and  energy  balance  in  a  tubular  flow  reactor 

Four  experimental  data  points  during  the  production  of 
nitroglycerine  in  a  tubular  flow  reactor  as  shown  in  Fig.  1, 
have  been  measured  by  Han  et  al.  (1988b).  These  include 
inflow  compositions  and  operating  conditions  and  are  listed  in 
Table  1.  We  can  evaluate  the  kinetic  parameters  of  the 
reaction  using  these  experimental  data. 


PROCESS  SAFETY  AND  ENVIRONMENTAL  PROTECTION  86  (2008)  37-47 


39 


J) 

1  F, 

[L() 

F  +dF 

/  \  r 

,  /  ^ 

system.  The  subscript;  of  Eq.  (10)  can  be  replaced  by  reactant. 
We  can  rearrangement  to  obtain 


Fig.  1  -  The  tubular  flow  reactor  of  producing 
nitroglycerine. 


If  we  consider  this  reaction  system  in  a  tubular  flow  reactor 
as  Fig.  1,  the  mass  balance  can  be  expressed  as 

Fg,o  dxG  =  (— rG)  dV  (5) 

where  FG,o  represents  the  molar  flow  rate  of  glycerine  and  V  is 
the  volume  of  tubular  flow  reactor,  respectively.  Substituting 
Eq.  (4)  into  Eq.  (5),  we  can  express  the  equation  of  mass  balance 
in  terms  of  fractional  conversion  of  glycerine  as 


=  Ae-E/RTCnG+m-1(l  -  xG)n(M  -  3xG)m  (6) 

where  r  =  V/vQ  is  the  space  time  and  v0  is  the  volumetric  flow 
rate  of  reactants. 

The  energy  balance  in  a  tubular  flow  reactor  is 
[Accumulation  of  heat]  =  [Enthalpy  flow  in] 

-  [Enthalpy  flow  out]  (7) 

At  steady  state,  there  is  no  heat  accumulated  in  this 
reactor.  Thus,  we  can  express  the  above  term  as 


(rG)hG(T)  +  (rN)hN(T)  +  (rNG)hNG(T)  +  (rw)hw(T) 

dT 

+  (FgC  p)G  +  FnCp;n  +  FnGCp;nG  +  FwCp,w  +  FsCp^)  =  0 

(11) 

The  subscript  G,  N,  NG,  W  and  S  represent  reactants  C3H5(OH)3, 
HN03,  C3H5(0N02)3,  H20  and  H2S04,  respectively.  The  correla¬ 
tion  of  reaction  rate  for  above  components  from  Eq.  (1)  is 


1  1 

-rG  =  -  ^  rN  =  -rw  =  rNG  =  r 


(12) 


Combining  Eqs.  (11)  and  (12)  one  obtains 

-  rhG(T)  -  3rhN(T)  +  rhNG(T)  +  3rhw(T) 

dT 

+  jjy  (FGCp,G  +  FnCp,n  +  FNGCP;NG  +  FWCP)W  +  FsCp,s)  =  0  (13) 


or 


r(AH)  +  gy  (FGCp)G  +  FnCp;n  +  FnGCp5nG  +  FwCp5w  +  FsCp;s)  =  0 

(14) 

The  generation  heat  in  the  glycerine-nitric  acid  reaction 
system  includes  the  heat  of  reaction  and  heat  of  dilution 
(Kaye,  1978b),  which  can  be  expressed  as 


0  = 


EhiT)  f,- 


h 1 


J  V+dV 


(8) 


where  h  is  the  enthalpy  of  reactant  which  can  be  expressed  as 
function  of  temperature  and  F  is  the  molar  flow  rate  of  reac¬ 
tant,  respectively.  Subscript  j  represents  the  reactants  of 
C3H5(OH)3,  HN03,  H2S04,  C3H5(0N02)3  and  H20  in  the  reaction 
system,  respectively.  Eq.  (8)  is  divided  by  dV  and  for  dV  — ►  0  it 
can  be  expressed  as 


_d_ 

dV 


Ehm 


j= i 


=  0 


(9) 


or 


AH  =  AHr  +  AHd  (15) 

where  AHr  represents  the  heat  of  reaction  and  AHd  represents 
the  heat  of  dilution,  respectively.  The  heat  of  reaction  also  can 
be  expressed  as 

A Hr  =  AHr,Tref  +  ACP  (T  -  Tref)  (16) 

ACp  =  Cp5nG  +  3Cp5w  —  Cp)G  —  3Cp5n  (17) 

where  Tref  and  AHr)Tref  are  the  temperature  and  heat  of  reaction 
at  standard  state.  Table  2  shows  the  values  of  heat  capacities 
and  heat  of  formation  at  standard  state  for  each  component  in 
this  reaction.  Substituting  Eq.  (16)  into  Eq.  (15)  we  can  obtain 


j=l  L 


d hj(T)  d T 
dT  dV 


=  0 


AH  =  AHr  +  AHd  =  AHr,Tref  +  ACP(T  -  Tref)  +  AHd 
(10)  =  AH0  +  ACP(T  -  Tref) 


(18) 


where  AH0  =  AHr;Tjref  +  AHd.  The  molar  flow  rates  of  C3H5(OH)3, 
where  dF/dV  =  q  represents  the  reaction  rate  of  reactant  and  HN03,  C3H5(0N02)3,  H20  and  H2S04  can  be  correlated  as  the 
dhj(T)/dT  =  CPj  is  the  heat  capacity  of  reactant  in  the  reaction  conversion  of  C3H5(OH)3  with  the  following  terms: 


Table  1  -  The  experimental  data  of  the  glycerine-nitric  acid  reaction  system  in  a  tubular  flow  reactor  [Han  et  al.  (1988b)] 


Test  code 


Inflow  compositions 


Operating  conditions 


Glycerine  Mixed  acid  Inflow  Outflow  Reaction 

- 7 — ; — 77 - 77 - 77 - 77 - 77  temperature  T0  (K)  temperature  Tf  (K)  time  r  (s) 

C3H5(OH)3  (g)  HN03  (g)  H2S04  (g)  H20  (g)  Amount  (g)  y  [  y  M  ;  W 


a 

32.72 

160.48 

562.08 

77.44 

800 

283.15 

307.65 

1.8 

b 

33.60 

192.96 

516.17 

94.87 

804 

282.95 

303.05 

2.1 

c 

33.17 

124.96 

567.28 

107.76 

800 

282.15 

301.15 

6.6 

d 

22.80 

68.34 

445.14 

86.52 

600 

282.35 

299.35 

14.0 
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Fg  =  FGjo(1  -  xG) 

(19) 

Fn  =  Fn,o  -  3FG;oXG 

(20) 

Fng  =  Fng,o  +  FG)oXG 

(21) 

Fw  =  Fw,o  +  3FG)qXG 

(22) 

Fs  =  Fs,o 

(23) 

Then,  substituting  Eq.  (5)  and  Eqs.  (18)-(23)  into  Eq.  (14),  we 
can  obtain 


dT 


AH0+ACP(T-Tref) 


Fg,o  dxG 


(Fg,oCp,g  +  Fn,oCp,n  +  Fng,oCp,ng  +  Fw,oCp,w  +  Fs,oCp,s) 
+  Fg,oxg(-Cp,g  -  3CPiN  +  CP;NG  +  3CP;w) 

CG)  o  dxG 

(CG!oCp;G  +  Cn,oCp,n  +  Cng,oCp,ng  +  Cw,oCp,w  +  Cs,oCp,s) 
+  CG!c>xG(-Cp)G  -  3Cp;n  +  Cp,ng  +  3Cp,w) 


Let  us  define 


(24) 


Cv  =  CG)oCp)G  +  Cn,oCp,n  +  Cng,oCp,ng  +  Cw,oCp,w  +  Cs,oCp,s  (25) 


Finally,  substituting  Eq.  (27)  into  Eq.  (6)  and  integrating  it, 
we  can  obtain  an  equation  for  glycerine  and  nitric  acid 
reaction  system  as  Eq.  (28). 


cG,rm,Cv  rT<  dT 

T  A(-AH)  JTo  e-E/RT[l  +  (Cv/CG,0AH)(T-T0)]n  1  1 

x  [M  +  3(Cv/Cg,o  AH)(T  -  To)]m 

In  general,  the  relationship  between  the  fractional  con¬ 
version  of  reactant  and  the  reaction  temperature  in  an 
adiabatic  reactor  can  be  expressed  as  x  =  (T  -  T0)/(Tf  -  T0)  on 
condition  that  the  reactant  is  entirely  reacted  during  the 
reaction  (Westerterp  et  al.,  1995a).  But  in  this  case,  the  reactant 
is  not  entirely  reacted  during  the  reaction.  Accordingly,  we 
derive  the  relationship  between  fractional  conversion  and 
reaction  temperature  for  this  reaction  system  as  Eq.  (27), 
which  can  substitute  into  Eq.  (6)  to  obtain  Eq.  (28).  The  data 
used  in  the  calculation  for  Eq.  (28)  are  listed  in  Table  3. 

2.2.  Analysis  of  heat  of  dilution  of  mixed  acid 


The  heat  produced  during  the  reaction  of  the  glycerine-nitric 
acid  system  is  not  only  due  to  the  heat  of  reaction,  but  also  due 
to  the  heat  of  dilution  of  the  acids  by  the  water  produced  in  the 
reaction.  The  heat  of  dilution  is  large  enough  to  influence  the 
reaction  temperature.  Therefore,  it  cannot  be  ignored  in 
calculating  the  total  enthalpy  (AF1).  The  heat  of  dilution  can  be 
determined  from: 

atj  WmfHdf  -  Wm0Hd0  ,  , 

AHd  = - Mng -  (29) 


Then,  substituting  Eq.  (25)  into  Eq.  (24)  and  integrating 
we  can  obtain 


L 


dT 


lTo  (AH0/ACP)  +  (T-Tref)  l  (Cv/CGi0ACp)+xG 


dxG 


Hdo  =  Hd,Tref  +  ACPd(T0  -  Tref) 


Hdf  =  Hd  Jref  +  ACPd(Tf  -  Tref) 


(30) 

(31) 


If  it  is  assumed  that  the  heat  capacities  and  heat  of  reaction 
are  constant,  then,  the  fractional  conversion  of  glycerine  can 
be  expressed  as  a  function  of  reaction  temperature.  Integrat¬ 
ing  results  in 


Cv(T-To) 
G  CGj0  AH 


(27) 


where  Wm0  and  Wmf  are,  respectively,  the  inflow  and  outflow 
masses  of  mixed  acid  in  reactor.  Hd0  and  Hdf  are  the  relative 
enthalpies  of  inflow  and  outflow  for  mixed  acid,  and  MNG 
represents  the  number  of  moles  of  nitroglycerine  after  the 
system  has  reacted.  Hd jTref  and  ACpd  are  relative  enthalpy  and 
specific  heat  of  mixed  acid  at  standard  state  (Tref=25  °C).  The 
values  of  relative  enthalpy  and  specific  heat  of  mixed  acid 


Table  3 

-  The  operating  conditions  of  glycerine-nitric  acid  reaction  system  in  a  tubular  flow  reactor 

Test 

code 

Reaction 
time  (t)  (s) 

Inflow  concentration  (mol/1) 

Outflow  concentration 

M  =  (Cn,o)/  (Cg,o) 

Cv  (kJ/lK) 

C3Hs(OH)3  (CGj0) 

HN03  (Cn, o) 

C3H5(OH)3  (CG,f) 

hno3  (CN,f) 

a 

1.8 

0.6258 

4.5274 

0.0063 

2.6689 

7.2347 

2.7029 

b 

2.1 

0.6163 

5.2209 

0.0062 

3.3904 

8.4711 

2.7176 

c 

6.6 

0.6346 

3.5265 

0.0063 

1.6417 

5.5570 

2.8182 

d 

14.0 

0.5948 

2.6297 

0.0059 

0.8632 

4.4213 

2.8738 
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Table  4  -  The  relative  enthalpy  and  specific  heat  of  mixed  acid  for  various  compositions  at  standard  state  (Tref  =  298.15  K) 
[Mckinley  and  Brown  (1942)] 


H2°  (p  , 

HNO3  ,p  X 

Relative  enthalpy  (kj/kg) 

Specific  heat  (kj/kg  K) 

hno3  +  H2SO4  +  h2o  1  w; 

HN03+H2S04 

30/100 

10/100 

-246.648 

2.039 

20/100 

-239.530 

2.094 

30/100 

-229.029 

2.207 

40/100 

-219.407 

2.290 

50/100 

-204.606 

2.315 

20/100 

10/100 

-213.617 

1.884 

20/100 

-205.485 

1.926 

30/100 

-198.079 

2.010 

40/100 

-187.858 

2.064 

50/100 

-173.995 

2.123 

15/100 

10/100 

-181.322 

1.830 

20/100 

-173.187 

1.872 

30/100 

-167.141 

1.901 

40/100 

-159.010 

1.943 

50/100 

-146.704 

2.010 

10/100 

10/100 

-134.122 

1.729 

20/100 

-127.348 

1.792 

30/100 

-123.282 

1.813 

40/100 

-119.731 

1.884 

50/100 

-113.371 

1.926 

0/100 

10/100 

-30.670 

1.453 

20/100 

-26.483 

1.537 

30/100 

-24.389 

1.578 

40/100 

-22.296 

1.620 

50/100 

-19.574 

1.645 

The  enthalpies  of  pure  H2S04,  HN03 

and  H20  are  assumed  as  zero. 

solution  for  various  compositions  were  obtained  from  Mckin- 

2.3. 

Determination  of  kinetic  parameters 

ley  and  Brown  (1942)  and  are  listed  in  Table  4.  We  can  express 
the  relative  enthalpy,  HdTref  and  the  specific  heat,  A Cpd  as 
functions  of  the  weight  percentage  of  water  in  mixed  acid 
Rw  and  the  weight  percentage  of  nitric  acid  in  anhydrous 
mixed  acid  RN  by  means  of  least  square  mathematical  curve 
fitting  approach  and  obtain  the  following  equations  as: 

HdJref  =;  -29.39  -  1276. 08RW  +  1689. 96R^  +  0.92RN 

+  296.24RwRn  +  42.74R2  (32) 

ACPd  =  1.4286  +  2.0398RW  +  0.5496RN  (33) 

We  can  use  Eqs.  (32)  and  (33)  to  calculate  the  values  of  the 
relative  enthalpy  and  the  specific  heat  in  the  ranges  of 
0  <  Rw  <  0.3  and  0.1  <  RN  <  0.5.  The  heat  of  dilution  for  four 
experiments  in  a  tubular  flow  reactor  are  calculated  and  listed 
in  Table  5.  It  is  found  that  the  experiment  ‘a’  has  a  larger  heat 
of  dilution  than  other  experiments,  which  leads  to  higher 
outflow  temperature  of  reactant  and  shorter  time  of  reaction 
in  this  experiment. 


The  kinetic  parameters  n ,  m,  A  and  E  of  the  glycerine-nitric 
acid  reaction  in  Eq.  (4)  can  be  determined  from  the  experi¬ 
mental  data  and  operating  conditions  presented  in  Tables  1 
and  3  according  to  Eq.  (28).  The  heat  of  dilution  calculated  and 
presented  in  Table  5  is  used  to  reduce  the  number  of 
unknowns  in  Eq.  (28)  to  4.  The  results  from  the  four 
experiments  are  used  to  determine  the  four  unknown 
parameters  with  the  integral  equations  being  solved  by  trial 
and  error  using  Mathematica  4.1  software  (Wolfram  Research, 
Inc.).  The  values  of  the  four  kinetic  parameters  n,  m,  A  and  E 
have  been  determined  to  be  0.9350, 1.117, 1.630  x  1021  (i1-n“m/ 
mol1_n-m  s)  and  122.0  (kj/gmol),  respectively. 

The  rate  of  glycerine  nitration  to  produce  nitroglycerine  is 
determined  by  the  concentrations  of  glycerine  and  nitric  acid, 
and  the  reaction  temperature.  The  temperature  of  the  reaction 
is  increased  by  both  the  heat  of  reaction  and  the  heat  of 
dilution  of  the  mixed  acid,  since  water  is  produced  during  the 
nitration.  Sometimes,  the  heat  of  dilution  exceeds  the  heat  of 
reaction.  The  heat  of  dilution  varies  with  the  reaction 
conditions,  which  thus  affect  the  rate  of  reaction.  In  this 


Table  5  -  The  dilution  heat  of  four  experiments  in  a  tubular  flow  reactor 

Test 

code 

Inflow  compositions  of  mixed  acid 

Outflow  compositions  of  mixed  acid 

AHd  (kJ/molNG) 

hno3  (g)  H2S04  (g) 

h2o  (g) 

Rw  (%) 

Rn  (%) 

hno3  (g) 

h2so4  (g) 

h2o  (g) 

Rw  (%) 

Rn  (%) 

a 

160.48  562.08 

77.44 

9.68 

22.25 

94.33 

562.08 

96.78 

12.85 

14.37 

-64.6104 

b 

192.96  516.17 

94.87 

11.8 

27.25 

125.03 

516.17 

114.73 

15.18 

19.50 

-47.5252 

c 

124.96  567.28 

107.76 

13.47 

18.05 

57.50 

567.28 

127.37 

16.93 

9.20 

-43.2737 

d 

68.34  445.14 

86.52 

14.42 

13.31 

21.97 

445.14 

100.00 

17.63 

4.70 

-42.0429 
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Fig.  2  -  The  relationships  between  reaction  temperature 
and  reaction  time  at  various  reaction  conditions  for  the 
glycerine-nitric  acid  system. 


x (sec) 


Fig.  3  -  The  relationships  between  fractional  conversion 
and  reaction  time  at  various  reaction  conditions  for  the 
glycerine-nitric  acid  system. 


investigation,  the  heat  of  dilution  was  only  one  determinant  of 
the  overall  heat.  The  relationships  between  both  reaction 
temperature  and  fractional  conversion,  and  reaction  time 
obtained  from  the  four  experiments  can  be  calculated  using 
the  values  of  reaction  kinetic  parameters  obtained  above. 
These  are  plotted  in  Figs.  2  and  3,  respectively.  Fig.  2  indicates 
that  reaction  temperature  varies  more  quickly  when  the  heat 
of  dilution  is  larger,  in  all  four  experiments  performed  in  a 
tubular  flow  reactor.  The  variation  in  fractional  conversion  is 
similar  to  that  in  Fig.  3.  The  numerically  simulated  reaction 


times  are  consistent  with  the  known  experimental  data 
obtained  from  four  experiments,  which  indicates  that  the 
kinetic  mode  of  derivation  in  this  investigation  is  reasonable 
for  a  glycerine-nitric  acid  reaction  system. 


3.  Safe  operating  conditions  of  nitroglycerine 
manufacture  in  Biazzi  continuous  process 

Although  nitroglycerine  has  been  produced  and  used  by  many 
countries  for  a  long  time,  the  production  process  is  still  very 
dangerous.  Numerous  thermal  runaway  incidents  and  explo¬ 
sion  disasters  have  occurred.  Unfortunately,  its  reactive 
hazard  has  not  been  clearly  identified.  Therefore,  operators 
must  carefully  maintain  the  entire  reaction  system  under 
proper  operating  conditions  to  prevent  any  runaway  of  the 
reaction  system  during  production.  The  Biazzi  continuous 
process  is  one  of  the  most  common  methods  for  nitroglycerine 
manufacture.  Two  factories  use  the  Biazzi  continuous  process 
to  produce  nitroglycerine  in  Taiwan.  Fig.  4  shows  a  general 
sketch  of  the  production  process  (Biazzi  NG  Plant,  1973).  The 
stable  reaction  criteria  and  safe  operating  conditions  in  this 
continuous  stirred  tank  reactor  (CSTR)  during  the  Biazzi 
continuous  process  are  determined  below. 

3.1.  Derivation  of  balance  equations 

In  the  first  place,  we  will  derive  the  mass  and  energy  balance 
equations  in  the  CSTR  for  glycerine-nitric  acid  reaction 
system.  The  transient  mass  balance  equation  for  the  produc¬ 
tion  of  nitroglycerine  in  the  CSTR  system  can  be  expressed  as 
follows: 

=  u0CG,0  -  u0CGf  -  V(-rG)  (34) 

where  —rG  is  the  reaction  rate  of  glycerine  which  can  be 
expressed  as  negative  sign  for  its  consumption.  The  u0  and 
V  are  the  flow  rate  of  reacted  solution  and  volume  of  reactant 
respectively.  CG)o  and  CG,f  are  inflow  and  outflow  concentra¬ 
tions  of  glycerine  in  this  CSTR. 

At  steady  state,  i.e.  dCc/dt  =  0,  dividing  Eq.  (34)  by  u0  and 
defining  the  resident  time  r  =  V/u0,  we  can  express  this 
equation  in  the  CSTR  as 

0  =  CG,0  -  CG,f  -  r(— rG)  (35) 

Substituting  Eqs.  (2)  and  (4)  into  Eq.  (35)  we  can  obtain  the  mass 
balance  in  the  CSTR  at  steady  state  as 

xG  =  rAe-(E/RTf)CG7"-1(l  -  xG)"(M  -  3xG)m  (36) 


Fig.  4  -  The  general  sketch  out  of  the  Biazzi  apparatus  (according  to  Biazzi  NG  Plant,  1973). 
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where  Tf  is  the  reaction  temperature  in  the  CSTR.  The  correla¬ 
tion  of  reaction  rate  for  C3H5(OH)3,  HN03,  C3H5(0N02)3  and  H20 
is  listed  in  Eq.  (12),  hence  the  outflow  concentration  of 
C3H5(OH)3,  HN03,  C3H5(0N02)3,  H20  and  H2S04  can  be 


expressed  as 

CG,f  =  Cg,o  -  rr  (37) 

CN,f  =  Cn,o  -  3rr  (38) 

Cngi  =  ^ng,o  +  rr  (39) 

Cw,f  =  Cw,o  +  3rr  (40) 

Cs,f  =  Cs,o  (41) 


where  A Hrjref  =  HNGjref  +  3Hw,Tref  -  HG,Tref  -  3HN)Tref  and  ACP  = 
Cp,ng  +  3CP)W  —  CP)G  —  3Cp>n. 

Substituting  Eqs.  (45)-(47)  and  (25)  into  Eq.  (44)  to  obtain 

CV(T0  -  Tf)  -  (2  AHr  +  AHd)rr  -  hS(Tf  -  Ta)/u0  =  0  (48) 

Substituting  Eqs.  (4),  (12)  and  (47)  into  Eq.  (48),  we  can  express 
the  energy  balance  in  the  CSTR  at  steady  state  as 

CV(T o  -  Tf)  -  [2  AHr,Tref  +  2  ACP(Tf  -  Tref) 

+  AHd]TAe-(E/RTf)C^+m(l  -  xG)n(M  -  3xG)m 

-hS(T  f-Ta)/u0  =  0  (49) 

To  analyze  the  equations  of  energy  and  mass  balance  in  a 
more  simple  form,  the  dimensionless  group  can  be  used  as 
following  definitions 


The  energy  balance  in  the  CSTR  is 
[rate  of  accumulation  of  heat  within  CSTR] 

=  [rate  of  heat  flow  into  CSTR] 

-  [rate  of  heat  flow  out  of  CSTR] 

+  [heat  generation  rate  by  chemical  reaction] 

-  [heat  removed  from  ambient  of  CSTR] 


¥>  =  *A  C£+m  1  exp  (—  ^) 

(42) 


At  steady  state,  there  is  no  heat  accumulated  in  the  CSTR. 
Thus,  the  above  term  can  be  expressed  as 


°°  =  RT2  (T°  “  Ta) 


(51) 


(52) 


(53) 


0  =  Uo(Cg,oHg,t0  +  Cn,oHn,t0  +  Cng,oHng,t0  +  Cw,oHw,t0  +  Cs,oHs,t0)  £ 

-  uo(CG,fHG,Tf  +  CNjfHNjTf  +  CNG)fHNG,Tf  +  CWjfHw,Tf  +  CSjfHs,Tf)  ^ref  =  Rp2  (^ref  _  ^a) 
+  V(-AHr-AHd)r-hS(Tf-Ta)  (43) 


where  Hg,t0  ,  HN;To ,  Hng,t0  ,  Hw,t0  and  HSj0  are  the  inflow  enthal¬ 
pies  of  components  C3H5(OH)3,  HN03,  C3H5(0N02)3,  H20  and 
H2S04,  respectively.  HGjf,HNjf,HNGjf,HWj{  and  HSj{  are  the 
outflow  enthalpies  of  the  above  five  components,  respectively. 
h  represents  the  convective  heat  transfer  coefficient  of  the 
cooling  medium,  S  denotes  the  external  surface  area  of  cooling 
coils,  Ta  is  the  inlet  temperature  of  cooling  medium.  Substitut¬ 
ing  Eqs.  (37)-(41)  into  Eq.  (43)  and  rearrangement  to  obtain 

Cg,o(Hg,T0  -  Hc,Tf)  +  Cn,o(Hn,T0  “  ^N,Tf)  +  CnG,o(HnG,T0  -  ^NG/Tf) 

+  Cw,o(Hw,T0  —  ffw,Tf )  +  Cs,o(Hs,T0  —  Hs,Tf) 

-  (HNG,Tf  +  3Hw,Tf  —  HG,Tf  —  3HN,Tf)rr 

+  (—A Hr  -  AHd)rr  -  hS(Tf  -  Ta)/u0  =  0  (44) 


-(2AHrJref  +  AHd)TAexp(-l/j8)Cgg" 
Cv/32(  E/R) 


(55) 


01  = 


-2ACPTAexp(-l/£)C£+m 


-G,0 


(56) 


hS 

UqCv 


(57) 


From  the  above  definitions  of  Eqs.  (50)-(57),  Eqs.  (36)  and 
(49)  can  be  simplified  into  the  following  equations 


Assuming  the  heat  capacity  of  each  component  is  setting  a 
reference  temperature  at  standard  state  as  Tref,  thus  we  can 
express  the  inflow  and  outflow  enthalpies  of  each  component 
as 

HjJo=CPJ(T0-Tref)  (45) 

HjJf  =CPJ(Tf-Tref)  (46) 

The  heat  of  reaction  is  expressed  as 

AHr  =  5>j,Tref  +  £(HjlTf  -  HjjJ  =  AHr,Trrf  +  ACp(Tf  -  Tref) 

(47) 


xG  =  pexp(j^-j)  (1  -  xg)"(M  -  3xG)m  (58) 

(4>o  +  fad  -  01e»ref)  exp  (^-j)  (1  -  xg)"(M  -  3xG)m 

=  (s  + 1)0  -  90  (59) 

Eqs.  (58)  and  (59)  are  the  dimensionless  mass  and  energy 
balance  equations  in  the  CSTR.  The  left  hand  side  and  right 
hand  side  terms  in  Eq.  (59)  represent  the  heat  generation  and 
heat  removed  rate,  respectively,  which  can  be  expressed  as 

Qg  =  (4>0  +  <£l<?  -  <Mref )  exp  (^-j)  (1  -  xG)n(M  -  3xG)m  (60) 


44 


PROCESS  SAFETY  AND  ENVIRONMENTAL  PROTECTION  86  (2008)  37-47 


and 

Qr  =  (e  +  1)0  -  6>o  (61) 

Thus,  Eq.  (59)  can  also  be  written  as 

Qg  =  Or  (62) 

In  general,  we  can  use  average  values  of  molar  heat 
capacity  Cp  and  reaction  heat  AHr  to  derive  the  energy  balance 
equation  for  a  CSTR.  We  can  take  either,  Cp  averaged  over  the 
feed  composition  in  the  temperature  range  inlet  to  outlet 
conditions  and  AHr  at  outlet  temperature,  or  Cp  averaged  over 
the  product  composition  in  the  same  temperature  range  and 
AHr  at  inlet  temperature  (Westerterp  et  al.,  1995b).  But  in  this 
research,  we  derive  the  energy  balance  equation  for  a  CSTR 
using  the  molar  heat  capacity  of  each  reactant  and  the  heat  of 
reaction  varies  with  the  reaction  temperature,  which  is  close 
to  the  actual  condition  of  the  reaction  system  compared  with 
above  mode  of  derivation.  Although  this  energy  balance 
equation  is  a  more  complicated  form,  we  use  the  dimension¬ 
less  groups  to  reduce  the  number  of  parameters.  This 
dimensionless  energy  balance  equation  can  aid  the  analysis 
of  safe  operating  conditions  for  the  glycerine-nitric  acid 
system  in  the  CSTR. 

3.2.  Verification  of  reaction  kinetic  parameters 

The  reaction  kinetic  parameters  are  estimated  from  the 
experimental  results  for  a  tubular  flow  reactor  using  the 
previously  described  numerical  analysis.  The  kinetic  para¬ 
meters  n,  m,  A  and  E  were  0.9350,  1.117,  9.780  x  1022  (ll  n  m/ 
rnol1-n-m  min)  and  122.0  (kj/g  mol),  respectively.  Tables  6  and 


Fig.  5  -  The  correlation  of  ln(-rG)  -  n  ln(CG>f)  and  ln(CN>f)  for 
glycerine-nitric  acid  reaction  system  in  a  CSTR  of  Biazzi 
continuous  process. 


7  present  the  operating  conditions  of  the  glycerine-nitric  acid 
reaction  system  in  a  CSTR  using  the  Biazzi  continuous  process 
(Klassen  and  Humphrys,  1953;  Biazzi  NG  Plant,  1973).  They  can 
be  used  to  verify  the  values  of  the  kinetic  parameters 
determined  above.  Fig.  5  plots  ln(-rG)  -  n  In  CG,f  versus  In  CN>f 
under  various  operating  conditions.  This  plot  yields  the  best 
linear  correlation  by  curve  fitting.  The  slope  and  point  of 
intersection  on  the  vertical  axis  of  the  line  equals  m  and  In  k, 
respectively,  where  k  =  A  e-E/RTf  is  the  kinetic  constant  of  the 
reaction  at  Tf  =  288.15  K.  The  values  of  m  and  In  k  are  1.117  and 
2.0122,  respectively,  which  correspond  to  the  results  of  the 
analysis  in  the  preceeding  section.  Accordingly,  the  kinetic 
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Table  7 

-  The  operating  conditions  of  glycerine-nitric  acid  reaction  system  in  Biazzi  continuous  process  (II) 

RF 

Rn 

Weight 
flow  rate  of 
inflow  (kg/min) 

Volumetric 
flow  rate  of 
inflow  (1/min) 

Inflow 

concentrations 

(mol/1) 

C3H5(OH)3  hno3 

Outflow 

concentrations 

(mol/1) 

c3h5(oh)3  hno3 

Resident  Reaction 

time  rate  (mol/lmin) 

(min) 

4.8 

49% 

6.148 

4.6265 

2.4904 

8.5535 

0.0290 

1.1696 

7.5650 

0.3254 

50% 

6.148 

4.6472 

2.4793 

8.6894 

0.0250 

1.3265 

7.5315 

0.3259 

51% 

6.148 

4.6678 

2.4683 

8.8240 

0.0219 

1.4847 

7.4982 

0.3263 

52% 

6.148 

4.6884 

2.4575 

8.9574 

0.0194 

1.6432 

7.4652 

0.3266 

4.9 

49% 

6.254 

4.7054 

2.4486 

8.5854 

0.0253 

1.3154 

7.4383 

0.3258 

50% 

6.254 

4.7264 

2.4377 

8.7216 

0.0221 

1.4747 

7.4051 

0.3262 

51% 

6.254 

4.7475 

2.4269 

8.8566 

0.0195 

1.6345 

7.3723 

0.3265 

52% 

6.254 

4.7686 

2.4162 

8.9904 

0.0175 

1.7943 

7.3397 

0.3268 

5.0 

49% 

6.360 

4.7842 

2.4083 

8.6163 

0.0223 

1.4585 

7.3157 

0.3261 

50% 

6.360 

4.8057 

2.3975 

8.7528 

0.0197 

1.6195 

7.2830 

0.3265 

51% 

6.360 

4.8272 

2.3868 

8.8881 

0.0176 

1.7805 

7.2506 

0.3268 

52% 

6.360 

4.8487 

2.3763 

9.0222 

0.0159 

1.9412 

7.2184 

0.3270 

5.1 

49% 

6.466 

4.8631 

2.3692 

8.6461 

0.0200 

1.5985 

7.1971 

0.3264 

50% 

6.466 

4.8850 

2.3586 

8.7830 

0.0179 

1.7608 

7.1648 

0.3267 

51% 

6.466 

4.9069 

2.3481 

8.9186 

0.0161 

1.9227 

7.1328 

0.3269 

52% 

6.466 

4.9288 

2.3376 

9.0531 

0.0146 

2.0841 

7.1011 

0.3271 

5.2 

49% 

6.572 

4.9419 

2.3314 

8.6750 

0.0182 

1.7353 

7.0823 

0.3266 

50% 

6.572 

4.9643 

2.3209 

8.8122 

0.0163 

1.8984 

7.0504 

0.3269 

51% 

6.572 

4.9866 

2.3105 

8.9481 

0.0148 

2.0610 

7.0188 

0.3271 

52% 

6.572 

5.0090 

2.3002 

9.0829 

0.0135 

2.2229 

6.9875 

0.3273 

parameters  of  the  glycerine-nitric  acid  reaction  system  can  be 
confirmed. 

3.3.  Calculation  of  heat  of  dilution 

The  heat  of  dilution  in  the  CSTR  can  be  assumed  to  be 

alt  WmfHdf  -  Wm0Hd0  .  . 

AHd= - Mng -  (63) 

where  Wm0  and  Wmf  are  weight  flow  rate  of  inflow  and  outflow 
of  mixed  acid  in  reactor,  respectively.  The  MNg  represents  the 
mole  producing  rate  of  nitroglycerine  after  the  system  is 


reacted.  The  Hd0  and  Hdf  are  inflow  and  outflow  relative 
enthalpies  of  mixed  acid,  respectively,  which  have  been 
described  by  Eqs.  (30)-(33).  We  have  derived  the  mass  balance 
equations  at  steady  state  in  the  CSTR  and  defined  as  Eq.  (36). 
The  fractional  conversion  of  glycerine  xG  at  various  reaction 
temperature  Tf  can  be  calculated  using  Eq.  (36)  by  means  of 
trial  and  error  technique.  Thus,  we  can  obtain  the  values  of 
dilution  heat  AHd  at  various  inflow  compositions  (RF,  RN), 
reaction  temperatures  (Tf)  and  inflow  temperatures  (T0)  using 
Eq.  (63)  in  the  CSTR  of  Biazzi  continuous  process.  Fig.  6  shows 
the  variations  of  AHd  at  various  Tf,  RN  and  T0  with  RF  =  4.8-5. 2 
for  glycerine-nitric  acid  reaction  system  in  the  CSTR  of  Biazzi 
continuous  process.  The  values  of  AHd  increase  along  the 


Fig.  6  -  The  variations  of  AHd  at  various  Tf,  RN  and  T0  with  RF  =  4.8-5. 2  for  glycerine-nitric  acid  reaction  system  in  a  CSTR  of 
Biazzi  continuous  process. 
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increasing  values  of  Tf,  but  decrease  along  the  increasing 
values  of  T0.  The  variations  of  RN  do  not  affect  the  value  of 
AHd.  Therefore,  the  dilution  heat  of  reaction  is  larger  when  the 
temperature  of  outflow  is  lower  or  the  temperature  of  inflow  is 
higher.  However,  the  variations  of  RF  affect  the  value  of  AHd, 
the  dilution  heat  of  reaction  is  larger  when  the  value  of  RF  is 
larger. 

3.4.  Analysis  of  optimum  operating  conditions  of  cooling 
system 

The  rate  of  chemical  reaction  is  very  sensitive  to  the 
temperature  of  an  exothermic  reaction  system.  When  the 
rate  of  generation  of  heat  in  the  reactive  system  exceeds 
the  rate  of  removal  of  heat  from  the  ambient  cooling  medium, 
the  whole  system  very  easily  accumulates  energy  and  its 
temperature  increases.  Then,  the  reactive  system  becomes 
unbalanced  and  a  fast  reaction  is  triggered.  As  soon  as  the 
temperature  reaches  the  critical  point,  the  reactive  system 
enters  an  unbalance  situation  or  even  explodes.  Therefore,  the 
heat  removal  rate  of  the  ambient  cooling  medium  must  be 
controlled  to  maintain  a  stable  temperature  of  the  reactive 
system. 

The  mass  and  energy  balance  equations  which  have  been 
derived  in  dimensionless  form  for  a  glycerine-nitric  acid 
reaction  system  in  the  CSTR  of  the  Biazzi  continuous  process 
are  shown  in  Eqs.  (58)  and  (59).  Tables  6  and  7  present  the 
operating  conditions  of  the  glycerine-nitric  acid  reaction 
system  of  the  Biazzi  continuous  process  in  the  CSTR;  these  can 
be  applied  to  analyze  the  optimum  operating  conditions  and 
its  cooling  system. 

Fig.  7  shows  how  the  fractional  conversion  of  glycerine,  xG 
is  affected  by  the  variations  of  0  and  p  for  RF  =  4.9  and  RN  =  51%. 
It  can  be  seen  that  the  fractional  conversion  of  glycerine 
increases  with  increasing  values  of  6  and  p,  this  means  that  a 
higher  temperature  of  the  reaction  system  or  a  higher 
temperature  of  the  cooling  system  will  lead  to  higher 
fractional  conversion. 

Fig.  8  shows  how  values  of  s  which  are  affected  by 
variations  in  0o  and  p  for  RN  =  49%  and  RF  =  4.8-5. 2  when  the 
temperature  of  reaction  system  is  set  at  15  °C.  e  is  the 
dimensionless  group  which  reflects  the  heat  transfer  capacity 
of  the  reactor.  It  can  be  seen  that  e  increases  with  increasing 


Fig.  7  -  Relationships  of  0-p-xG  for  glycerine-nitric  acid 
reaction  system  in  a  CSTR  of  Biazzi  continuous  process. 


Fig.  8  -  Relationships  of  0o-p-s  at  RN  =  49%  for  glycerine- 
nitric  acid  reaction  system  in  a  CSTR  of  Biazzi  continuous 
process. 


Fig.  9  -  Relationships  between  Ta  and  T0  at  RN  =  49%  and 
RF  =  4.8-5. 2  for  glycerine-nitric  acid  reaction  system  in  a 
CSTR  of  Biazzi  continuous  process. 


values  of  0o  and  p,  which  means  that  the  use  of  a  higher 
reactant  inlet  temperature  or  a  higher  temperature  in  the 
cooling  medium  requires  higher  values  of  h,  the  overall  heat 
transfer  coefficient  or  S,  the  surface  area  of  the  cooling  coils. 
The  relationships  of  90-p-s  at  other  values  of  RN  are  similar,  but 
the  value  of  s  is  larger  when  the  value  of  RN  is  larger. 

There  are  two  factories  which  use  the  Biazzi  continuous 
process  to  produce  nitroglycerine  in  Taiwan.  Their  operating 
conditions  are  considered  from  the  viewpoints  of  both  safety 
and  economy.  The  inflow  rate  of  glycerine  and  mixed  acid  are 
controlled  at  1.06  and  5.194  kg/min,  respectively.  The  weight 
ratio  of  mixed  acid  to  glycerine  is  4.9  and  the  weight 
percentage  of  nitric  acid  in  mixed  acid  is  51%.  The  temperature 
of  inflow  compositions  is  kept  within  18-20  °C.  The  tempera¬ 
ture  of  the  nitrator  is  maintained  at  15  °C  by  cooling  coils  fed 
with  brine.  The  cooling  surface  and  heat  transfer  coefficient  of 
coils  are  1.7  m2  and  47.104  kJ/K  m2  min  (Biazzi  NG  Plant,  1973), 
respectively.  The  inlet  temperature  of  the  cooling  medium  is 
adjusted  within  -2  to  -5  °C  in  the  nitrator.  We  can  use  the 
above-mentioned  system  parameters  to  predict  the  relation¬ 
ship  between  the  inlet  temperature  of  cooling  system,  Ta  and 
the  inlet  temperature  of  the  reactants  entering  the  reactor,  T0 
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with  compositions  within  RF  =  4.8-5. 2  and  RN  =  49%  from  Fig.  8. 
The  result  is  shown  in  Fig.  9.  The  relationships  of  Ta  and  T0  at 
other  values  of  RN  are  similar,  but  the  value  of  Ta  is  larger  when 
the  value  of  RN  is  smaller  under  the  same  value  of  T0.  From  the 
prediction  shown  in  Fig.  9,  we  find  the  inlet  temperature  of 
cooling  system  must  be  kept  within  -4.55  to  -5.05  °C  when  the 
inlet  temperature  of  the  reactants  T0  is  18-20  °C.  The  result  of 
this  simulation  is  close  to  the  actual  operating  conditions  of 
the  reaction  system. 


4.  Conclusions 

From  the  above  calculations  and  analyses,  this  research 

obtains  the  following  conclusions: 

1.  The  reaction  rate  of  glycerine-nitric  acid  to  produce 
nitroglycerine  in  the  CSTR  of  Biazzi  continuous  process  can 
be  expressed  as  -rG  =  9.78  x  1022  e_(14674  044/Td  C^35^1/7, 
which  has  been  verified  from  several  actual  operating 
conditions. 

2.  The  heat  of  dilution  of  reaction  is  varied  by  operating 
conditions  for  glycerine-nitric  acid  reaction  system  in  the 
CSTR  of  Biazzi  continuous  process.  The  heat  of  dilution  is 
larger  when  the  temperature  of  outflow  is  lower  or  the 
temperature  of  inflow  is  higher. 

3.  The  relationship  of  60-/3-b  at  various  values  of  RN  and  RF  is 
shown  in  Fig.  8,  which  can  be  used  to  evaluate  the  heat 
transfer  coefficient  of  reactor  which  is  an  important 
parameter  for  the  reactor  design  of  glycerine-nitric  acid 
reaction  to  produce  nitroglycerine. 

4.  The  relationship  between  inlet  temperature  of  cooling 
system  Ta  and  temperature  of  inflow  reactants  T0  at  various 
values  of  RN  and  RF  is  shown  in  Fig.  9.  These  relationships 
can  be  used  for  safe  operation  in  the  present  factories  of 
Taiwan. 
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